Hidden Topological Order in 3 Na (F = 1) Bose-Einstein Condensates 
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We show the existence of a new hidden topological order in 23 Na (F=l) Bose-Einstein condensates 
(BEC) with antiferromagnetic interactions. Occurrence of this order is due to the confinement of 
hedgehogs in spin ordered BEC where a spin Josephson effect takes place. However, a topological 
long range order is also argued to coexist with a short range spin correlation, as a result of topological 
order from disorder. 



23 Na(F = 1) BEC with antiferromagnetic interactions 
turn out to have many extremely fascinating properties 



One of them is the topology of the order param- 
eter space. It was pointed out in previous works that 
23 Na BEC in optical traps are in Quantum Spin Nematic 
States, with the order parameter space being [S 1 x S 2 ]/ Z 2 
P,P[. In this article we will report a new hidden long 
range order of pure topological nature and explore the 
possibility of having spin correlated states based on a 
topological consideration. To study spin correlated BEC 
of spin-1 interacting atoms, it is important to properly 
choose a set of collective coordinates. These variables are 
to completely characterize the nonlinear spin dynamics 
in the whole parameter space and to yield a simplest rep- 
resentation for two-body spin dependent scatterings. For 
this reason we employ a vector n living on a unit sphere 
for the description of the nonlinear spin dynamics in the 
presence of antiferromagnetic interactions. The low en- 
ergy spin dynamics in the system can be mapped into an 
o(3) nonlinear sigma model(NLerM) Rpf, 



£=-(cyi) 2 ,n 2 = l, 

/ = (16 7 r) 1 / 2 ( /0 Aa 3 ) 1/6 ,Aa 



a 2 - a 



(1) 



ao,2 are the scattering lengths in F = 0, 2 channels and 
p is the density. We also have introduced dimension- 
less length and time: r — > r/p 1 ' 3 , it — ► r/vsp 1 ' 3 , and 

v s = JlQitTi Aap/M 2 . Derivatives <9 M are defined as 
(d T ,d x ,d y ,d z ). A reduction from [S 1 x S 2 ]/Z 2 to S 2 is 
possible when Z 2 fields discussed in M are effectively 
frozen, or Z 2 strings are gapped. Moreover, the local 
coupling between the U(l) phase and spin order involves 
higher derivative terms and becomes important only at 
a rather high frequency Kp 2 ' 3 /2M. 

Following Eq.l, the energy of the system consists of 
two parts: a) the potential energy H 2 p(Vn) 2 /2M, which 
is the energy cost in the presence of a slow variation of n, 
b) and the zero point kinetic (rotation) energy Io(dtn) 2 /2, 
where Iq = (h Aap/M)^ 1 can be considered as the ef- 
fective inertial of an individual atom. This inertial origi- 
nates from two-body scatterings and is inversely propor- 
tional to scattering lengths. Zero point rotations tend 
to disrupt the order of n between different atoms. / is 



a measure of the amplitude of quantum fluctuations in 
23 Na BEC. 

In weakly interacting limit / is much less than unity, 
the symmetry is broken in the ground state and the BEC 
have a spin order, with n(r, t) = no. Just as the Joseph- 
son effect occurs in a superfluid where a phase long range 
order is established, a spin Josephson effect occurs in spin 
ordered BEC. A spinor condensate in a weak external 
magnetic field Be z exhibits the following time depen- 
dence: 

*(*) = -=[exp^+t)|l,l > +fflq>(v»_t)|l,-l >], (2) 

and p± = pi ± u>Bj ^b = 9^bB (g is the Lande factor 
of atoms). Notice that in the presence of a spin stiffness, 
a well defined relative phase between two components of 
the condensates is developed. To observe the spin co- 
herence, one introduces a reference condensate which is 
at a chemical potential p 2 but experiences no magnetic 
fields. The temporal dependence of the interference be- 
tween these two condensates would be 



5p(t) <x cos ((/ii — p 2 )t) cos (u>Bt). 



(3) 



At u>b — 0, Eq. 3 reflects to a usual AC Josephson effect 
and when p\ — p 2 , one observes an AC spin Josephson 
effect due to a spin order in BEC, a sinusoidal-time de- 
pendence with frequency determined by wg instead of 
difference in chemical potentials. The alternative way to 
observe the spin coherence without introducing a refer- 
ence condensate is to transform |1, 1 > component into 
|1, — 1 > after the magnetic field is switched on for time 
t. This results in an interference between the induced 
1 1,-1 > component and the original |1, — 1 > compo- 
nent. The condensate density in the overlapped region 
exhibits a sinusoidal dependence on the waiting time t 
with a frequency 2lub- 

Quantum fluctuations of the nematic order can be 
estimated in a lowest order approximation. Introduc- 
ing <5n = n — no and assuming the fluctuations are 
weak, we obtain, < (<5n 2 ) >oc /. As / increases, spin 
wave excitations start to interact strongly and Eq.2 be- 
comes invalid. The renormalization group (RG) equa- 
tion of / is determined by the interactions between col- 
lective modes. In 3 + 1 dimension, the RG equation 



takes a form df/dl = (3(f), with the /3-function given as 
(3(f) = -2/(/ c - /). Within the frame work of NLcrM, 
f c = 8ir 2 in d = 3. / characterizes spin correlations in 
the ground state of 23 Na BEC. 

At a high density limit f > f c JlQ], zero point kinetic 
energy dominates and the spin stiffness is renormalized 
to zero at a long wave length limit. Especially, in an 
extremely quantum disordered phase, the potential en- 
ergy at interatomic scale Et — h p 2 ! 3 J2M is negligible 
compared with zero point rotation energy E = h /2Iq 
of an individual atom, n of each atom fluctuates inde- 
pendently and spins of atoms only correlate at an inter- 
atomic distance. The excitation spectrum has a gap of 
order E . As the two-body scattering gets weaker, the 
inertial of each individual atom Iq is higher and the mo- 
tion of n on the unit sphere becomes slower. Thus, n 
of different atoms starts to correlate at a finite length £, 
which is a function of the parameter /, much longer than 
the interatomic distance. Under the current situation, 
£ = /0 _1 / 3 / c /[/ — fd- And the state is rotation invariant. 
The temporal correlation is characterized by T£, the time 
needed for a spin wave to propagate over a distance £ 
with a velocity v s . At a low density limit / < f c , the 
potential energy dominates and the spin stiffness flows 
to a finite value under the renormalization group trans- 
formation. All spin-1 Bosons rotate as a rigid body with 
an effective inertial I(N) = NIq. The energy gap in 
the excitation spectrum vanishes E gap ex h /2NIq as the 
number of Bosons N approaches infinity, signifying the 
existence of gapless Goldstone mode. 

For 23 Na BEC confined in highly anisotropic traps, 
zero point motions of collective variables are more promi- 
nent. Let us approximate highly anisotropic traps as ID 
systems. The /3-function RG equation in this case can 
be shown as (3(f) — l/27r/ 2 , which always flows into a 
strong coupling point, i.e. a rotation invariant state. And 
the correlation length is given as £ = p" 1 exp(47r//) (p 
is the linear density). However, the length of traps has 
to be longer than £ in order for this spin liquid to be 
observed in experiments. 

To study the topological long range order, we intro- 
duce a connection field to characterize all configurations 
of n [Ql . The field strength is a gauge invariant Pontrya- 
gin density F M „ = e^e^n^^n,,^^. In (3 + 1)D, 
Foi, i = X,y,x corresponds to an electric field and 



e ijk ¥ 



Hfe represents three components of magnetic 



field; in 2d, Foj, i = x,y is an electric field with two 
components and the magnetic field has only z compo- 
nent H z = F xy . Particularly, we will be interested in the 
fluctuations of a topological charge defined as 



' a 



— IdSH- 

4tt J s 



(4) 



where the integral is carried over the boundary of volume 
V and e„ is a unit vector normal to the boundary. We 
are going to show that at large V, t limit, 



lim < (C m (t) - C m (0)Y 



>=0 



(5) 



in both spin ordered and spin weakly disordered BEC, 
as results of a topological long range order. However, in 
strongly disordered limit, 



9 tV 

lim < (C m (t) - C m (0)) >= — 

t,V^oo T m 



(6) 



where r m is finite depending on tunneling rate of mag- 
netic monopoles(see below). Averages in Eqs.5,6 are 
taken over the many-spin wavefunction ^({n}). In 2D, 
we define C m = (47r) _1 f dxdyU z , the total number of 
Skyrmions living on the 2D sheet and Eqs.5,6 hold in 
spin ordered and spin disordered BEC respectively (V is 
replaced by an area S) . The topological long range order 
occurs when hedgehogs are confined, either due to the 
condensation of atoms or due to an order from disorder 
mechanism. 

Let us start with 2D. First, we notice that in space 
x = (r, r), H,, = e w,y F M j, is the solution of the Gauss 
equation 



9 r; H r; (x) = VQ m <5(x-x m ) 



(T) 



in the presence 
Euclidean space. 



of space-time monopoles {x™ 1 } in (2 + l)d 



(a) 

mimmti 

(b) 



Fig.l. Spin configurations at a) t = +oo and b) t = — oo 
boundaries of a monopole in (2 + l)d. 

The space-time monopoles represent instantons, which 
essentially connect a trivial vacuum to a Skyrmion con- 
figuration, as seen in Fig.l. In Euclidean space, a space- 
time monopole located at x — (tq,Tq) terminates a 
Skyrmion centered at ro at time to. Each monopole 
causes a change in topological charge C m by one unit; the 
topological density C m — 1 is conserved in the absence 
of space-time monopoles in (2 + 1)D. Indeed, following 
Eqs.4,7, 



dC m (r) _ ^ 



rS(T 



(8) 



where the surface contribution has been neglected since 
we will be interested in C m per unit volume. 



The probability of finding deconfined space-time 
monopoles depends on the energy of Skyrmions with 
C m = 1 with respect to C m = configuration. In po- 
lar coordinates (r, <fi) , a static Skyrmion is a configura- 
tion with n(p,cf>) — (sinO(r) cos </>, s'm6(r) sin0, cos(r)); 
0(r) varies from at r = to i at r > (, with £ 
an arbitrary parameter. For spin ordered BEC where 
spin Josephson effects should be observed, the energy 
of the Skyrmion is proportional to 8ttH p/2M and is 
scale invariant. For a Skyrmion of given £ |12| , the con- 
nection field is concentrated in a region of the size £, 
H z = r _1 sin 6(r) 89 (r)/dr, and gets screened at a length 
scale larger than £. In spin ordered BEC, a Skyrmion 
is nondegenerate with respect to a trivial vacuum. As 
we will see, this leads to the confinement of space-time 
monopoles. 

To illustrate the point of confinement, we consider C m 
as a function of time at some discrete unit and find the 
following binary representation for monopoles. 



a) ...00000000000000000.. 






y 



b) ...00000001111111111. 



c) ...00000111111100000... 



Fig. 2. In a binary representation, monopole-like instan- 
tons are represented by kinks living on a string, with C m , 
the topological charge as the order parameter. 

Each bit shown in Fig. 2 labeled as 0, or 1 represents 
the topological charge read out at certain moment; a 
1-bit corresponds to a Skyrmion while a 0-bit is for a 
C m = configuration. The a) string represents a trivial 
vacuum where C' m remains zero at all time. In the second 
binary string b), a 1-string stands for any scale invari- 
ant Skyrmion created at the time to defined by the do- 
main wall separating 0-string and 1-string. So a domain 
wall in our binary representation represents a monopole 
Q m = 1, terminating a Skyrmion configuration (C m = 1) 
at the interface and changing C m by one unit. A domain 
wall-anti domain wall pair in the c) line corresponds to a 
monopole- antimonopole pair with Q m = ±1; a 1-string 
here, or a Skyrmion, is terminated at a monopole at one 
end t\ and at an antimonopole at the other end £2- 

Since the energy of a Skyrmion in spin ordered BEC 
is higher than that of a trivial configuration, the energy 
of a 1-bit is positive with respect to a 0-bit. The energy 
of a domain wall in Eq.l which is the number of 1-bits 
in the structure, defines the action of the monopole and 
is proportional to the length of the 1-string. The ac- 
tion of an isolated space-time monopole, or creation of a 
Skyrmion therefore is proportional to L t fi p/2M, being 
infinity. (L t is the perimeter along a temporal direction.) 



The action of having monopole-anti-monopole pairs is 
proportional to (£1 — £2), the time interval between cre- 
ation and annihilation of Skyrmions, thus monopoles are 
confined. For this reason, Eq.8 vanishes and a topological 
long range order prevails. 

The conclusion arrived so far can be extended to 3D 
straightforwardly. The change of the topological density 
is due to the quantum nucleation of monopoles instead 
of Skyrmions. In the binary representation in Fig. 2, a 1- 
bit stands for a static monopole and a 0-bit for a trivial 
configuration. A domain-wall represents a termination 
of a static monopole at certain time. Since the energy 
of a monopole is proportional to the system size L in 
spin ordered BEC, the monopole is nondegenerate with 
respect to the trivial vacuum. Following the same argu- 
ment carried out in 2D, the action to have a monopole 
nucleated is proportional to LL t h p/2M, and is infin- 
ity. This result has a profound impact on the monopole 
confinement. The interaction between a monopole and 
an antimonopole has to be linear in terms of distance 
between them. The action of a monopole pair nucleated 
at t\ and annihilated at £2 with largest spatial separation 
\r\ — r<i\ is proportional to |j"i — 7*2 1 (ti — £2) and monopoles 
are confined because of spin stiffness. Therefore the right 
hand side of Eq.8 has to be zero because of the confine- 
ment, and topological charge C m is a conserved quantity. 
We once again arrive at Eq.5 in a spin ordered BEC and 
topological charge is conserved. 

The situation in disordered limit is more delicate and 
depends on dimensionalities. In spin disordered BEC, 
the spin stiffness is renormalized to zero and spin fluc- 
tuations are gapped. In the absence of spin stiffness, 
the energy of a Skyrmion is determined by its interac- 
tion with spin fluctuations. Upon integration of spin 
wave excitations, based on a general consideration of the 
gauge invariance and the parity invariance, we conclude 
Eq.l should be reduced to £ S (F M „) = F^ u F^/2g + ...; 
g~ 1 (A s ) is a function of the spin gap A s measured in 

units of (E T E ) 1/2 = H 2 (Aap 5 / 3 ) 1/2 /2M. In a leading 
approximation, g~ l (x) is a logarithmic function of x, i.e. 
\tlx~ 1 in 3-D and is a linear function of x in 2D p3[ . 
We will analysis the topological order based on C s . The 
energy of C m -configuration is aiVjCl^ if quantum tun- 
neling is neglected, with a(L) a function of the system 
size. 

In 2D, with induced interactions, Skyrmion energy is 
no longer scale invariant and is minimized at £ = 00. The 
connection field of a Skyrmion spreads over the whole 2D 
sheet. The energy of a Skyrmion scales as L~ 2 and van- 
ishes as the system size L goes to infinity. This implies 
that skyrmion configurations become degenerate with 
a trivial vacuum and monopoles are deconfined. This 
mechanism of deconfinement resembles the liberation of 
fractionalized quasiparticles in one dimension polymers 
pT| . In one dimensional polyacetylene, the ground state 



has twofold degeneracy because of the Peierls instability 
and domain wall solitons become free excitations. 

When instantons are liberated, the energy of a con- 
figuration C m is ill-defined. A more serious consider- 
ation beyond the mean field approach outlined above 
involves the evaluation of monopole-like instantons ac- 
tion. In (2 + \)d the partition function of a monopole 
configuration {x m }, is V^ n P"/rdexp( — Sin). And 
P = L 2 L t exp(— w/g) is the quantum tunneling ampli- 
tude; S in = g^ 1 Ysm'^m ^QmQm' /|x m ' - x m | represents 
the Coulomb interaction between space-time monopoles. 
The above result suggests that Skyrmions always con- 
dense in the spin disordered BEC and topological charge 
C m is not a conserved quantity. It is important to real- 
ize that in the absence of Zi strings, ±1 Skyrmions can 
be physically and homotopically distinguished because 
of the coherence of quantum spin nematic BEC, unlike 
the situation in a classical nematic liquid crystal. By 
moving a Skyrmion around a 7r-disclination, the conden- 
sate acquires a 7r-phase with respect to the original one, 
which can manifest itself in a Josephson type of effect. 
One can also distinguish ±1 Skyrmions by looking at the 
local connection field. In a positive Skyrmion configu- 
ration, a spin-i collective excitation, which carries half 
charge with respect of connection fields, experiences a 
connection field of an opposite sign compared to that of 
a negative Skyrmion Q . 

So we are able to show that the change of C m has 
a short range temporal correlation because of instanton 
effects, 



< 



dCm(T) dC m (0) 



Ot 



Ot 



>=cxp( ) 



(9) 



And T m is proportional to £ 2 T£ expfjr/g). Immediately, 
one recognizes that Eq.9 simply indicates a random walk 
of C m as a function of time in spin disordered 2D BEC. 
Therefore , we conclude Eq.6 holds. 

But in 3D weakly disordered BEC where stiffness has 
vanished, a static monopole still carries a finite energy 
because of spin fluctuations induced interactions. In 313, 
as illustrated in £ S (F M „), spin fluctuations discriminate 
topologically different configurations. Particularly, the 
fluctuations are strongest in C m = configuration so 
that the energy is lowest in the topological trivial con- 
figuration. The situation differs from 2D spin disordered 
limit where the energy of topological nontrivial configu- 
ration with C m = 1 is the same as that of C m — one. 
For a monopole, each spherical shell of radius R can be 
viewed to be a Skyrmion squeezed into a 2D sheet of 
size R. The smaller shells, or textures of a finite size, 
dominate the energy of a monopole. Spin fluctuations 
interacting with the singularity lift a degeneracy between 
the monopole configuration and a trivial vacuum in the 
absence of spin stiffness. It also appears to us that this is 
independent of the form of C s (F fi „) we introduced. The 



energy of a monopole here is inversely proportional to g, 

1 /2 

i.e. (E Et) lnA.T 1 . The action of having a monopole, 
which leads to a change in C m by one unit, is proportional 
to L t and is infinite for this reason. 

So the monopoles remain confined even after spin cor- 
relation becomes short ranged, with the action to have a 
pair of monopoles being proportional to (ti —t<z). Only 
pair production which conserves C m is allowed. The 
topological long range order thus coexists with a short 
range spin correlation in 3D. This can be considered as 
a case of order from disorder phenomena. And, the con- 
finement of monopoles in this case is driven purely by the 
spin fluctuations, instead of spin stiffness in spin ordered 
BEC. However, by increasing the spin disorder further, 
the degeneracy between C m — 1 and C m — could be 
established, which signifies deconfinement of monopoles 
and breakdown of the long range topological order. 

It is also possible to demonstrate the topological or- 
der/disorder in terms of space-time correlations of con- 
nection fields or Wilson loop integrals. The topologi- 
cal order is vital for the quantum number fractionaliza- 
tion of excitations. In 3D spin disordered BEC, where 
a topological long range order is established, spin-i ex- 
citations are elementary ones |1J]. On the other hand, 
when hedgehogs are liberated or condense, only spin-I 
excitations exist in the spectrum. Finally, quantum hid- 
den orders in spin liquids of strongly correlated electrons 
are recently reviewed by Wen |l5| . Topological orders in 
spin triplet superconducting liquids are also investigated 
by Dernier et al. @]. I would like to thank ASPEN 
center of Physics and Amsterdam summer workshop on 
"Flux, charge, topology and statistics" for their hospi- 
talities. Discussions with A. Abanov, Duncan Haldane, 
N. Read, T. Senthil, O. Starykh and Paul Wiegmann are 
greatly acknowledged. 
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